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I.  INTRODUCTION 


1. 1  This  volume  of  the  comprehensive  report  presents  the  specific 
computer  simulation  that  has  been  devised  and  tested  during  the  contract 
year, 

1.2  In  Section  II  the  conventional  computer  flow  charts  are  presented. 
In  Section  III  the  symbols  used  in  the  flow  charts  are  defined  and  a  de¬ 
scription  of  the  flow  charts  is  given.  This  description  takes  the  reader 
step-by-step  through  the  simulation  explaining  the  purpose  of  each  com¬ 
puter  operation. 

1.3  In  Section  IV  the  literal  computer  program  for  the  simulation  is 
presented.  In  order  to  achieve  generality  in  application,  two  versions 
of  the  program,  one  in  ALGOL  and  the  other  in  FORTRAN,  have  been  pre¬ 
pared  and  are  presented  one  after  the  other. 
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FLOW  CHARTS  FOR  COMPUTER  SIMULATION  OF 
ARTILLERY-FIRED  GB  SHELLS 


2. 1  This  section  presents  the  seven  flow  charts  which  provide  a  con¬ 
ventional  representation  of  the  computer  simulation. 

2.2  With  the  advent  of  the  new  computer  languages  such  as  ALGOL 
and  FORTRAN,  certain  simplifications-have  become  available  in  computer 
programing  which  tend  to  cause  minor  differences  between  the  flow  charts 
and  the  programs  based  on  them.  For  example,  the  FOR. instruction  in 
ALGOL  supplants  the  conventional  loop  program  illustrated  in  the  flow  charts. 
The  instruction:  FOR  N  =  (1,  1,  12);  BEGIN. ..  .takes  the  place  of  the 

loop:  SET  N  =  1;  PERFORM  OPERATION  ON  VARIABLE  N;  Is  N  =  12?;  IF 
N  =  12,  EXIT;  IF  N  /  12,  ADD  1  TO  N,  PERFORM  OPERATION  ON  VARIABLE N. ... 

2.3  When  several  loops  are  contained  one  inside  the  other,  the  flow 
charts  tend  to  clarify  the  order  of  operations  which  can  become  obscure 
when  embodied  in  the  computer  programs. 

2.4  Fiow  Chart  1  gives  an  over-all  view  of  the  simulation.  The  suc¬ 
ceeding  charts  present:  (a)  targeting  and  delivery  of  the  GB  shells,  (b); de¬ 
termination  of  the  fragmentation  effect,  (c)  determination  of  the  splash 
effect,  (d)  determination  of  the  GB  vapor  effect.  (This  las.  appears  on 
three  charts.  On  the  first  chart  of  the  series  a  geometric  determination  of 
whether  a  target  is  in  the  path  of  the  GB  vapor  cloud  is  made.  On  the 
second  chart  the  determination  of  breathing  parameters,  response  delay 
time,  breathing  phase  at  impact  time,  masking  time,  breathhold  time,  and 
gasp  time  is  made.  On  the  third,chart  the  inhaled  vapor  dose  is  calculated.) 
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2.  TARGETING  AND  DELIVERY 
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III.  DESCRIPTION  OF  THE  FLOW  CHARTS 


3. 1  This  section  presents  a  list  of  symbols  used  in  the  flow  charts 
and  a  step-by-step  description  of  each  operation.  Each  explanation  is 
keyed  to  a  corresponding  block  on  the  flow  charts  by  the  use  of  block 
numbers . 

3.2  Although  Flow  Chart  1  (the  Master  Flow  Chart)  provides  for  de¬ 
termination  of  time  dependent  impingement  (downwind  impingement),  there 
is  no  flow  chart  or  program  at  present  to  carry  out  this  determination.  This 
feature  of  the  model  has  not  yet  been  developed,  but  is  expected  to  be 
added  early  in  the  coming  contract  year.  The  modular  nature  of  the  model 
makes  this  addition  quite  simple. 

3.3  Since  this  model  is ,  by  design,  a  Monte  Carlo  simulation,  it  is 
of  interest  to  summarize  the  variables  whose  values  are  determined  by, 
sampling  from  probability  distributions.  In  the  current  form  of  the  simula¬ 
tion,  the  following  probability  distributions  are  used. 

a.  Target  element  location  error  (assumed  symmetrical 
in  x  and  y  !/). 

b.  Battery-volley  centroid  x-coordinate  error. 

c.  Battery-volley  centroid  y-coordinate  error. 


The  assumed  symmetry  in  x  and  y  makes  the.detec'  on  process  correspond 
to  aerial  photographic  reconnaissance.  If,  for  example,  artillery  ob¬ 
server  reconnaissance  were  used,  there  would  be  different  distributions 
for  x  and  for  y. 
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d.  Individual- piece  x-coordinate  error. 

e.  Individual-piece  y-coordinate  error. 

f.  Response  delay  time. 

g.  Masking  time. 

h.  Duration  of  breathhold. 

1.  Tidal  volume. 

3.4  With  further  development  it  will  be  possible  to  sample  values  of 
additional  quantities  such  as  source  strength,  stability  parameter,  and 
wind  speed. 

3.5  In  ad.j.r:-.*,  to  the  sampling  of  probability  distributions  just  re¬ 
ferred  to.,  there  ere  currently  three  other  random  determinations: 

a.  tit  the  start  there  is  a.  random  selection  of  three 
target  elements  frormthe  array  of.target  elements 
(see  Block  T—  1) . 

b.  In  the  fragmentation  routine  a  target  may  be  determined 
to  be  in  a  zone  of  potential  fragmentation  casualty. 

When  this  'is  the  case  a  random  number  is  compared  with 
the  probability  of  casualty  for  the  particular  zone.  If 
the  random  number  is  less  than  or  equal  to  the  probability 
the  target  is  a  casualty  {see  Block  B— 16). 

c.  In  the  breathing  parameter  subroutine  the  target's 
phase  in  the  breathing  cycle  is  selected  at  random. 

3.6  Currently,  the  chemical  dose-response  relation  is  not  used  in 
the  program  which  has  been  written.  The  final  output  of  a  computer  run 
is  the  total  chemical  dose  received  by  each  target  element.  It  is  a  rela¬ 
tively  simple  matter  to  extend  the  program  so  that  a  Monte  Carlo  deter¬ 
mination  of  degree  of  casualty  can  be  made  following  the  determination 
of  the  chemical  dose. 

V 

NOTATION  FOR  PLOW  DIAGRAM 

3.7  The; notation  utilized  in  the  flow  diagram  is  defined  under  the 
following  classifications: 

1;  Coordinate  Notation 

2.  Targeting  and:  Delivery  Notation, 


3.  Fragmentation  Notation 


4.  Splash  Notation 

5.  Breathing  Parameter  Notation 

6.  Dosage  Notation 


Coordinate  Notation 

(Xj,  Yj) 

<V  yn> 

(x1,  y') 


(X,  y) 

(x",  y") 

(x^,  yiw) 

(xw  vw) 
'  n  'yn' 

(x1,  y\) 


Location  of  impact  of  the  itn  shell 

Location  of  the  n  target 

Location  of  one  of  the  targets  observed 

Estimated  location  of  one  of  the  targets  observed 

.Estimated  centroid  of  target  observed; 
location  at  which  battery  fires 

Location  of  impact  centroid 

Errors  in  firing  of  individual  guns 

Location  of  impact  of  i**1  shell  in  wind  coordinates 

frVi 

Location  of  n  target  in  wind  coordinates 

Distance  downwind  and  crosswind  from  impact  to 
target. 


Targeting  and  Delivery  Notation 


Impact  index;  i  1.2..... i * 

Location  of  impact  of  shell  fired  from  left-most 
gun  if  there  were  no  individual  piece  firing  error. 

Target  index;  n  =  1 , , .  •  #  nMAX* 

The  largest  integer  obtained  by  dividing  the  number 

1000  by  nMAX‘ 

The  distance  between  aiming  points  of  each  gun. 

The  direction  of  the  wind  in  degrees  measured 
counterclockwise  from  the  x-axis. 


Fragmentation  Notation 


Distance  of  the  n*b  target  from  the  itb  impact. 

Cotangent  of  the  angle  between  the  x-axis  and 
the  line  between  the  impact  and  the  target. 
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A 

A  small  distance  from  the  impact  in  which  a  maximum 
casualty  occurs. 

^MAX 

Maximum  range  of  fragmentation. 

T 

Sector  In  which  fragmentation  occurs;  r  =  1,2, . 

/* 

c 

Denotes  varying  radii  as  measured  from  the 
impact  joints;  ^  1,2,.,., £j^AX(t) 

«rl 

Cotangent  of  the  angle  between  side  1  of  the  rtb 

sector  and  the  x-axis . 

°V2 

Cotangent  of  the  angle  between  side  2  of  the  rtb 
sector  and  the  x-axis. 

fT* 

The  radii  for  the  sector. 

ZT* 

Splash  Notation 

The  zone  formed  by  the  sector  and  the  radii, 

r7£andr7,£-l 

k 

Index  for  denoting  ellipses  of  different  sizes; 
k  —  1,2,..,,  k 

ak 

Distance  from  the  center  of  the  kth  ellipse  to  the 
edge  of  the  ellipse  along  the  major  axis. 

bk 

Distance  from  the  center  of  the  k*b  ellipse  to  the 
edge  of  the  ellipse  along  the  minor  axis. 

ck 

Distance  from  the  center  of  the  kl“  ellipse  to  each 
focus . 

*k 

Distance  from  the  center  of  the  k**1  ellipse  to 
the  impact  point. 

Pnik 

Sum  of  the  distances  of  the  target  from  each  of  the 
two  foci  for  the  ktb  ellipse. 

Breathina  Parameter  Notation 

J 

Gasp  index;  j  =  1 , 2 , . . . , j^AX* 

BHTj 

Time  for  jtb  breathhold  plus  exhale  time. 

GD 

Gasp  Duration  (duration  of  inhale  portion  of  gasp). 

q 

Indexes  an  interval  of  width  6  in  the  breathing 
cycle;  q  = 
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R 


Burst  radius  of  shell. 
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*A 

^m 

tm 

*w 

t-1 


*E' 

'Gi' 

tHl' 


TV 

a,b 


Increment  of  time  in  breathing  cycle  usedfor 
numerical  integration. 

Time  of  arrival  of  gas  cloud  relative  to  time  of 
impact. 

Time  of  departure  of  gas  cloud  relative  to  time 
of  impact. 

Masking  time  (in  increments  of  6). 

Time  at  which  masking  is  completed  relative  to 
time  of  impact. 

The  time  relative  to  the  beginning  of  a  man's  breathing 
cycle  at  impact  time. 

Response  delay  time  (increments  of  6). 

Time  relative  to  time  of  impact. 

Time  when  gas  cloud  arrives  relative  to  beginning 
of  the  man's  breathing  cycle. 

Time  at  which  the  jth  "gasp"  begins  relative  to 
beginning  of  the  man's  breathing  cycle 

Time  at  which  breathholding  first  begins  relative 
to  beginning  of  breathing  cycle. 

Time  at  which  masking  is  completed  relative  to 
beginning  of  the  man’s  breathing  cycle. 

The  time  relative  to  the  beginning  of  the  man's 
cycle  when  the  peak  of  the  breathing  cycle  is  obtained. 

The  length  of  the  breathing  cycle. 

Time  relative  to  beginning  of  the  breathing  cycle 
when  the  man  can  respond  to  the  impact. 

Tidal  volume  of  air  intake. 

Integers;  the  duration  of  inhale  is  2a6;  the  duration 
of  exhale  is  2b6;  tp'  =  2a6;  tj'  -  tp'  =  2b6. 

Location  for  accumulating  steady  breathing  dose. 
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q 

u 

a 

y 

rf(q) 

£v(q) 

X^x'y’) 


Location  for  accumulating  gasp  dose.. 

Intervals  in  breathing  cycle;  q  =  1,2,  .  ..tj'^. 

Wind  speed. 

Stability  parameter. 

Index  used  for  numerical  integration. 

Fraction  of  TV  inhaled  during  the  q**1  interval. 

th 

Volume  of  air  intake  during  the  q  interval. 

Concentration  of  cloud  at  time  t  when  target 
is  x'  and  y'  from  impact. 


DESCRIPTION  OF  FLOW  CHARTS-ACC-2,  20  JUNE  1961 


Master  Flow  Chart 

BLOCK  M-l.  Enter  Pilot  Model 

The  input  elements— the  constants  and  probability  distributions  necessary 
for  the  battle-field  model  considered— are  entered  into  the  computer.  • 

BLOCK  M-2.  Input  Target  Elements 

The  (x,  y)  coordinates  of  the  targets  are  entered  into  the  computer.  The 
number  of  target  elements  varies  from  3  to  100. 

BLOCK  M-3.  Select  Aiming  Point 

The  point  at  which  the  guns  are  to  be  fired  is  determined  as  ,a  consequence 
of  randomly  selecting  three  of  the  actual  target 
with  target  acquisition  are  considered)-. 

BLOCK  M-4 .  Compute  Impact  Points 

The  coordinates  of  each  of  the  impact  points  are  determined.  Two  sets  of 
coordinate  systems  are  used  througiiout  the simulation 

1.  a  regular  coordinate  system  and 

2.  a  wind  coordinate  system. 

'8LOCK.M-5 .  Js  Frag  Considered? 


elements  (errors  associated 


This  block  makes  consideration  of  fragmentation  effects  optional. 

BLOCK  M-6.  Compute  f’raq  Effects/Casualties 

This  block  determines  the  fragmentation  effects. 

BLOCK  M-7 .  Is  Splash  Considered? 

This  block  makes  consideration  of  splash  effects  optional. 

BLOCK  M-8.  Compute  Splash  Dose/Casualties 

This  block  determines  the  splash  effects. 

BLOCK  M-9.  Is  Time  Dependent  Vapor  Considered? 

This  block  makes  consideration  of  time  dependent  vapor  effects  optional. 

BLOCK  M-10.  Compute  Vapor  Dose/Casualties 

This  block  determines  the  results  of  the  time  dependent  vapor. 

BLOCK  M-ll.  Is  Time  Dependent  Impingement  Considered ? 

This  block  makes  consideration  of  time  dependent  impingement  effects 
optional. 

BLOCK  M-12.  Compute  Impingement  Dose/Casualties 

This  block  determines  the  results  of  the  time  dependent  impingement. 
BLOCK  M-13.  Cumulate  Chemical  Dose 

This  block  determines  the  total  chemical  dosage  obtained  from  the  splash 
effect,  the  time-dependent  vapor  effect,  and  the  time-dependent  impinges 
ment  effect. 

BLOCK  M-14.  Compute  Net  Chemical- Casualties 

This  block  determines  the  total  chemical  casualties  as  a  result  of  the 
chemical  dose  obtained  in  Block  M-13. 
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BLOCK  M-15. 


Is  There  to  Be  Another  Replication? 

This  block  tells  the  computer  the  number  of  times  to  perform  the  simulation. 

If  there  is  to  be  another  replication,  the  sequence  of  events  that  began 
in  Block  M-l  is  repeated. 

BLOCK  M-l 6.  Stop 

When  all  replications  for  the  specified  input  elements  are  completed,  the 
computer  stops. 

The  blocks  as  generally  described  in  the  Master  Flow  Chart  will  be  described 
in  more  detail  in  the  following  pages . 

Targeting  and  Delivery 

BLOCK  T-l.  Select  at  Random  Three  Targets 

From  the  n^  target  elements  entered  into  the  machine ,  throe  target  elements 
are  selected  at  random.  The  random  process  for  selecting  these  three  targets 
is  as  follows: 

1.  The  largest  integer,  r',  that  can  be  obtained  by  dividing  the  number 
1000  by  is  determined. 

2.  Using  the  interval  width  of  r\  intervals  are  constructed. 

3.  There  is  an  interval  which  corresponds  to  each  of  the  target  elements, 
l,2,...,rij (Any  remainder,  1000  -  (r’Hn^j^),  as  a  consequence 
of  dividing  nj^^  into  1000,  is  not  a  usable  ranaom  number.) 

4.  The  first  target  element  is  obtained  by  determining  the  correspondence 
between  the  interval  in  which  the  first,  random  number  falls  and  the 
input  target  elements .  If  the  randpm  number  selected  is  greater  than 
(r'MnMAX>,  another  random  number  is  chosen. 

5.  After  the  first  target  is  selected,  there  are  two  intervals  which  are  not 
usable,  the  interval  corresponding  to  the  first  target  selected  and  the 
interval  greater  than  (r')(njv|ftx)- 

6.  This  process  is  repeated  until  three  different  target  elements  have  been 
selected. 

BLOCK  T-2.  Calculate  the  Centroid 

The  centroid  of  the  three  targets  selected  at  random  is  calculated  by  the 
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following  equations: 


El  x<*/3  y=  zl  V3 4 

O'  =1  O'  =  1 

BLOCK  T-3.  Estimate  Coordinates  of  Randomly  Selected  Targets 

The  x  and  y  errors  involved  in  the  estimates  of  the  three  target  elements 
are  determined  from  a  probability  distribution  of  errors .  In  the  current 
version  of  the  model  the  variance  of  the  x  errors  is  identical  with  the 
variance  of  the  y  errors.  These  variances  can  be  different.  The  estimated 
coordinates  of  the  three  randomly  selected  targets  are  obtained  by  adding 
these  x  and  y  errors  to  the  actual  x  and  y  coordinates  of  the  three  randomly 
selected  target  elements ,  respectively. 

BLOCK  T-4.  Calculate  Estimated  Centroid 

The  estimated  centroid  is  the  point  at  which  the  center  of  the  volley  is 
aimed.  This  point  is  obtained  from  the  estimated  coordinates  obtained 
in  Block  T-3  by  the  following  two  equations: 

3  3 

*=  Z]V/3  y  =  2Z' y«/3 

a  -  1  a  =  L 

BLOCK  T-5.  Print  Targeting  Information 

The  following  targeting  information  is  , printed: 

1.  The  coordinates  of  each  of  the  randomly  selected  target  elements  - 
(determined  in  Block  T-l). 

2.  The  actual  centroid  of  the  three  randomly  selected  elements 
(determined  in  Block  T-2). 

3.  The  coordinates  of  the  three  estimated  target  elements  (de¬ 
termined  in  Block  T-3). 

4.  The  estimated  centroid  for  firing  (determined  in  Block  T-4). 

BLOCK  D-l .  Select  Impact  Centroid 

The  impact  centroid  is  determined  by  adding  the  x  and  y  c-jrors  in  battery 
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firing  to  the  coordinates  of  the  estimated  centroid  (calculated  in  Block  T-4). 

The  x  and  y  errors  are  obtained  from  a  probability  distribution  of  errors  in 
the  impact  centroid  from  the  estimated  centroid.  The  variance  of  the.x 
errors  is  not  equal  to  the  variance  of  the  y  errors.  This  block  takes  into 
account  the  artillery  battery  delivery  error. 

BLOCK  D-2.  Select  Impact  Errors  for  Each  Impact 

The  x  and  y  errors  for  each  impact  point  about  the  point  at  which  the  gun 
was  fired  are  obtained  from  a  probability  distribution  of  individual  piece 
firing  errors .  An  x  and  y  error  is  selected  for  each  of  the  i^j^  impact 
points . 

BLOCK  D-3.  Computer  Impact  Coordinates  in  Regular  (x.  vl  System 

Figure  1  illustrates  the  computation  of  the  impact  coordinates.  The  x-coordinate 
of  the  impact  point  of  the  shell  fired  from  the  left-most  gun  would  be  located 
a  distance  of  (i  v*ax“1)(C/2)  =  I*  from  the  x-coordinate  of  the  impact  centroid 
(determined  in  Block  D-l)  if  there  were  no  individual  piece  firing  error.  The 
remaining  x-coordinates  of  the  points  at  which  the  shells  would  strike  if  there 
were  no -individual  piece  firing  errors  are 

L  +  ?  ,  L+  2?,...,L+  Umax-I)?. 

By  adding  the  x- errors  as  obtained  in  Block  P-2  to  L  +  £,  L  +  2£, » . . , 

L  +  (ij^r-1)?,  the  actual  x-coordinates  of  the  individual  impact  points 
are  determined.  The  actual  y-coordinates  for  each  of  the  Individual  impact 
points  are  determined  by:  adding  the  y-errors,  as  obtained  in  Block  D-2, 
to.the  y-coordinates  of  the  impact  centroid,  as  obtained  in  Block  D-l. 

BLOCK  D-4  .  Compute  Wind  Coordinates  of  Impact  Points 

The  impact  coordinatesare  determined  when  the  regular  x  and  y  axes  are 
rotated  by  the  number  of  degrees  corresponding  to  the  wind  direction. 

The  following  equations  are  used: 

xi  =  xicos®  +  ViSinB 
y**  =  yjCosB  -  XjSinO 

where  x.  and  are  the  impact  coordinates  in  the  regular  coordinate 
system  (where  the  x  axis  is  parallel  to  FEBA  and  the  y  axis  is  perpen¬ 
dicular  to  FEBA),  xw  and  yw  are  the  coordinates  in  the  wind-coordinate 
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FIGURE  1.  COMPUTATION  OF  THE  IMPACT  COORDINATES 
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system,  and  0  is  the  wind  direction  in  degrees  with  respect  to  the  x  axis 
(the  xw  axis  is  parallel  to  wind  direction) . 


BLOCK  D-5.  Print.  Delivery  Information 

The  impact  centroid  and  the  impact  points  for  each  of  the  shells  are 

printed  both  in  the  regular  coordinate  system  and  also  in  the  wind  coordinate 

system. 

Fragmentation 

BLOCK  F-l.  Is  Frag  Considered? 

This  block  makes  consideration  of  fragmentation  effects  optional. 

BLOCK  F-2.  Set  i  =  1  and  n  =  1 

In  this  block  i  and  n  are  set  equal  to  one. 

BLOCK  F-3.  Is  y,  y  ? 

'1  •'n 


See  Figure  2.  The  coordinates  ,  of  the  impact  point  are  denoted  by  (xj.yj), 
and  the  coordinates  of  the  target  are  denoted  by  (xn,yn).  Whenever  the 
y  coordinate  of  the  impact  point,  yj,  is  great?  *  than  the  coordinate  of  the 
target,  yn,  the  target  is  not  hit. 

BLOCK  F-4.  Compute  Distance  of  Target  from  Impact  Point 

For  economy  in  computer  operation,  the  squares  of  distances  are  computed 
and  compared,  rather  than  the  actual  distances.  The  distance  squared 
of  the  target  from  the  impact  point  is  computed  by  the  following  equation: 

dni2  =  (xn  -  */  +  (yn  -  y/  .  , 

BLOCK  F-5.  Is  dni2  P  MAX  ? 


There  is  a  distance,  from  the  impact  point  beyond  which  fragmenta¬ 

tion  effects  are  negligible.  If  fragmentation  effects  are  negligible,  the 
computer  returns  to  Block  F-21.  If  the  target  is  within  the  severity 
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FIGURE  2.  FRAGMENTATION  PATTERN 


of  fragmentation  effects  is  computed  in  the  remaining  blocks. 

BLOCK  F-6.  Is  dni8  s  ? 

It  is  possible  that  the  distance  between  the  target  and  the  impact  point 
is  so  small,  that  a  direct  hit  is  obtained.  If  a  direct  hit  is  obtained, 
the  occurrence  of  such  a  hit  is  stored. 

BLOCK  F-7 .  Store  Maximum  Casualty 

If  a  direct  hit  has  occurred,  maximum  casualty  is  stored. 

BLOCK  F-8.  Is  yd  =  yn? 


In  this  block,  and  also  Block  F-9,  and  F-10,  attention  is  focused  oh 
determining  the  cotangent,  j3nj,  of  the  angle  between  thex  axis  and  the 
line  connecting  the  impact  point  and  the  target.  If  the  impact  point  ana 
the  target  are  the  same  distance  from  the  x-axis,  yj  =  y  ,  fsn^,  is  infinite. 

BLOCK  F-9.  Set  0ni  =  co 


Since  the  computer  is  not  equipped  to  handle  a  number  as  large  as  infinity, 
special  handling  must. take  place.  The  location  which  contains  /?ni  is  set 
equal  to  9.9  x  1040  to  indicate  that  0ni  is  equal  to  infinity, 

BLOCK  F-10,  Compute  /3ni 


Since  this  cotangent  is  less  than  infinity,  the  computation  for  /?ni  is 
obtained  by  using  the  following  equation: 


BLOCK  F- 1 1 .  Set  r-»  1 

T  denotes  the  sectoi  in  which  fragmentation  effects  occur.  (See  Figure  2.) 
With  T  ~  1,  a  determination  w'ill  be  made  to.  see  if  the  target  is  in  the 
ffrstsector.  Ir,  genera  1,  there  are  sectors  (r  =  1,2, . . .  .r^^x). 


BLOCK  F-12. 


Set  £  =  1 


The  Greek  letter  £  is  used  as  a  subscript  to  denote  varying  radii  as 
measured  from  the  impact  point.  In  general,  £  is  equal  to  1,2,...  *£maX(t). 
The  radii  for  the  sector  are  denoted  by  rTi,  rr2, . . .  <rT£iyi AX(r)‘ 

BLOCK  F-13.  Is  <yTi  £  Sni? 


As  shown  in  Figure  2,  aT\  is  the  cotangent  of  the  angle  between  side  one 
of  the  r™  sector  and  the  x  axis.  If  is  not  less  than  or  equal  to  f}^, 
the  target  is  not  in  the  sector.  In  this  case  the  machine  returns  to 
F-21  and  the  next  impact  is  considered,.  If  ctji  is  less  than  or  equal  to 
pn |,  a  check  is  made  in  Block  F-14  against  side  two  of  the  t**1  sector  . 

BLOCK  F-14.  .  Isj9ni£(yT2? 


The  notation  a^z  is  used  to  denote  the  cotangent  of  the  angle  between 
side  two  of  the  rth  sector  and  the  x  axis.  If  0ni  is  equal  to  or  less  than 
aT2,  the  target  is  In  the  r**1  sector. 

BLOCK  F-15.  Is  dni  2  :>  rT£  2  ? 


Since  £  is  initially  set  to  one,  a  no  answer  means  that  the  target  lies 
within  radius  one  of  sector  t.  This  region  is  denoted  zone  Zr  j.  If. the 
answer  is  yes,  £  is  increased  until  the  zone  is  determined. 

BLOCK  F-16.  Calculate  Fragmentation  Effects  for  Zone  ZT  ^ 

Corresponding  to  each  zone,  ZT^,  there  is  a  probability  that  a  frag¬ 
mentation  casualty  occurs.  A  casualty  occurs,  if  a  random  number  is 
equal  to  or  less  than  the  probability  of  a  casualty.  The  number  of 
casualties  for  each  target  element  is  accumulated  over  all  impacts. 

BLOCK  F-17 .  Is  £  equal  to  £  MAX(t)^ 


If  £  is  equal  to  £mAX(t)'  the  tar9®t  element  is  in  an  area  free  from  frag¬ 
mentation.  The  sequence  of  operations  in  F-21  follows  and  the  effects 
of  the  next  impact  are  considered.  If  £  is  not  equal  to  £maX{t)/  ‘hen 
the  target  element  may  be  within  a  radius  corresponding  to  a  greater  value 
of  £  and  the  next  larger  radius  is  considered. 
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BLOCK  F-18. 


Increase  g  by  one 


For  the  same  sector  of  r,  £  is  Increased  to  the  next  radius,  £  +  1.  The 
sequence  of  operations  in  Block  F-15  follows. 

RLOOK  F-10.  Tst  =  t . ? 

.  ~  '  MM 


If  T  is  equal  to  r  MAX,  then  the  target  element  is  outside  the  last  fragmen¬ 
tation  sector  and  fragmentation  effects  are  negligible.  If  t  is  not  equal 
to  tMAX  {that  is'  r  is  less  than  tMAX^'  th®  next  T  sector  is  considered. 

BLOCK  F-20.  increase  t  by  one 

Since  t  is  less  than  not  all  sectors  have  been  considered.  There¬ 

fore,  the  next  sector  is  considered  by  increasing  t  by  one.  After  increasing 
r,  the  sequence  of  operations  in  Block  F-12  follows. 

BLOCK  F-21.  Clear  dfii2  and  /3nj 


Since  the  next  impact  is  to  be  considered,  certain  storage  locations  must 
be  cleared. 

BLOCK  F-22.  Is  i=  ij^? 


If  i  is  equal  to  iw^,  all  impact  points  for  this  target  element  have  occurred 
and  a  print-out  follows.  If  i  is  not  equal  to  ij^  (that  is,  i  is  less  than 
iMAX^'  further  impact  points  are  to  be  considered. 

BLOCK  F-23 .  Increase  1  by  one 

The  next  Impact  Is  considered  by  increasing  i  by  one. 

BLOCK  F-24.  Print  Out 

The  results  of  the.  Impacts  on  the  target  elements  are  printed  out. 
Specifically,  the- number  of  times  in  which  the  target  element  is  a  casualty 
is  printed  out. 

BLOCKF-25.  Set  1  =  1 

of  the  impacts  on  the  next  target  element  are 
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considered,  again  beginning  with  the  first  impact. 
BLOCK  F-26.  Is  n  =  nj^? 


If  n  is  equal  to  njyjAX'  the  effects  of  all  impacts  on  each  of  the  target 
elements  have  been  evaluated  aha  printed;  therefore,  the  computer  enters 
the  routine  for  evaluating  the  effects  pf  splash.  If  n  is  ndt  equal  to 
nMAX  is'  ”  is  i®ss  t^ian  nMAx)'  the  effects  on  all  of  the  target 
elements  have  hot  beer,  evaluated. 

BLOCK  F-27.  Increase  n  by  One 

The  value  of  n  is  increased  by  one  so  that  the  fragmentation  effects  on 
the  next  target  element  may  be  evaluated. 

Splash 

BLOCK  S-l.  Is  Splash  Considered? 

This  block  makes  consideration  of  splash  effects  optional. 

BLOCK  S-2.  Set  i- 1,  n«  1,  k-k^ 


This  block  initializes  values  for  the  three  parameters  1,  n,  and  k.  The 
letter  k,  is  used  as  an  index  for  denoting  ellipses  of  different  sizes; 
as  k  increases,  the  size  of  the  ellipse  increases. 

Note  to  BLOCKS  S-3,  S-4,  S-5,  and  -6: 

In  Blocks  S-3,  S-4,  S-5,  and  S-6,  the  index  k  has 
only  the  value  kj^^  so  that  only  the  largest  splash 
region  is  under  consideration.  The  largesttsplash 
ellipse,  kj^jj>  is  approximated  by  a  circumscribing 
rectangle.  If  the  target  element  lies  inside  this 
rectangle,  further  tests  starting  with  BLOCK  S-7 
are  employed.  These  tests  determine  the  annular 
region  between  ellipses  in  which  the  target  element 
lies. 

BLOCK  S-3.  Is  yn  >  yt  +  Ak  +  bk? 


This  block  determines  whether  the  target  element  Is  too  far  downrange 
from- the  point  of 'Impact. 
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BLOCK  S-4.  Is  yn  <  yj  +  Xk  -  bk? 


This  block  determines  whether  the  target  element,  is  too  far  uprarigejfrom 
the  point  of  impact. 

BLOCK  S-5.  Is  xp  >  xi  +  ak? 


This  block  determines  whether  the  target  element  is  too  far  to  the  right 
of  the  point  of  Impact. 

BLOCK  S-6.  Is  xn  <  xi  -  a^? 


This  block  determines  whether  the  target  element  is  too  far  to  the  left 
of  the  point  of  impact. 

BLOCK  S-7.  Compute.  pnlk 


Refer  to  Figure  3.  The  symbol  Pnik  is  used  to  denote  the  sum  of  the 
distances  of  the  target  from  each  of  the  two- foci,  of  the  kt^  ellipse  and 
ck  is  used  to  denote  the.distance  from  the  center  of  the.k— \ellipse  to 
each-focus.  The  formula  for  pnlk.is.as  follows: 

Pnik  =1  {xn  -  xi)Z  +  *  n+  ck  -V2  +{k  ~  xi)2  +  ^  ~ck 54  xk>2  • 


BLOCK  S-8. 


Is  Pnik  *  Za.kMAX? 


If  p.j,  >  2akjvj^jj«  then  the..target.  element ds  not  within  the  ellipse 

and  tne  next  impact  point  will-,  be.  considered.  If  p^IL,  is  less/t^  " 


then  the  target- element  is  within  the  k 
blocks  determine  the  lowest  numbered  e 
element. 


a.i,  jo  iMBHiian  23],  .  , 

ellipse..  The  subsequent 
"pse  which  ;  includes  the  target 


BLOCK  S-V. 


Is  k  .=  1? 


If  k  -is  equal  to  one ,  the. target  element  is  within  the  first.ellipse  or  region 
one.  -Proceed  to  .Block  S.- 13  .  Ifk.is  not  equal.to  one,  proceed  -to  Block 
Sr  10, 
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Pntk  =  length  of  I  + 
length  of  II 


FIGURE  3.  SPLASH  PATTERN  SHOWING  KM.V  ELLIPSE 
AND  CIRCUMSCRIBED  RECTANGLE 
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BLOCK  S-10. 


Decrease  k  by  One  and  Compute  pn^ 


The  value  of  p  ^  for  the  next  smaller  ellipse  is  computed. 

BLOCK  S-ll.  I§  Pnik5,  2ak? 

iL 

If  p  ^  is  greater  than  2a^,  then  the  target  element  is  not  within  the  kin 
ellipse,  but  is  within  the  k  +  I**1  ellipse.  An  inquiry  is  made  about  the 
target  being  in  the  k4*1  ellipse  only  if  the  target  is  within  the  k  +  l**1 
ellipse.  If  p  is  less  than  28^,  then  the  target  element  is  within  the 
k*h  ellipse.  Therefore,  the  target  element  is  in  one  of  the  following 
regions:  k,  k  -  1, . . . ,  2,  1. 

BLOCK  S-12.  Determine  Splash  for  Region  k  +  1 

Since  it  has  been  ascertained  that  the  target  element  is  within  the  k  +  l**1 
region,  the  corresponding  amount  of  splash  is  determined.' 

BLOCK  S-13.  Determine  Splash  for  Region  1 

As  a  result  of  the  logic  in  Block  S-9,  it  is  determined  that -the  target 
element  is  within  region  one  and  the  corresponding  quantity  of  splash 
is  determined. 

BLOCK  S-14.  Set  k  =  k.„.v 

. _ MAX 


Since  the  location  of  the  target  element  with  respect  to  the  spiash  regions 
has  been  determined  for  this  impact,  the  value  of  k.is  reset  to  ^mAX  *9r 
use  in  determining  the  location  of  the  target  element  with  respect  to  the 
next  impact. 

BLOCK  Sr  15 Is.  i  =  ij^  ? 


If  i  is  equal  to  iMAX»  all  impacts  have  been  considered  with  reference  to  a 
specific  target  element.  If  i  is  not  equal  to  iw»v  (that  is,  i  is  less  than, 
Imax)'  th®n  hext  impact  is  to  "be  considereaT 

■BLOCK:S=16.  Increase  1  by- One 

The  next -impact  is  considered  by  increasing  i  by  one. 


BLOCK  S-17. 


Is  n  ~  nMAX? 


If  n  is  equal  t°  njyjAX'  tar<3st  elements  have  been  evaluated  for  all 
impacts  and  the  computer  will  proceed  to  the  evaluation  of  time-dependent 
vapor.  If  n  is  not  equal  to  (that  is,  n  is  less  than  nMAX)»  the  next 

target  element  is  considered. 

BLOCK  S-18.  Increase  n  by  one,  set  i  =  i 

In  order  to  ascertain  the  effects  of  the  splash  upon  additional  target 
elements  for  each  of  the  impacts,  n  is  increased  by  one  and  i  is  set 
equal  to  one . 


Time  Dependent  Vapor 

BLOCK  V- 1 .  Is  Time  Dependent  Vapor-Considered? 

This  block  makes  consideration  of  time  dependent  vapor  effects  optional. 


BLOCK  V-  2 .  Set  1  =  1:  n  =  1 


This  block  initializes  values  for  the  parameters  i  and  n. 


BLOCK  V- 3.  Calculate  x*  and  y’  distances  from  Impact  to  target 

(See  Figure  40  This  block  calculates  the  x’  distance  and  y‘  distance  (wind 
coordinates)  between  the  target  and  the  impact.  The  equations  used  for 
this  calculation  are  as  follows: 


„  w  w 

A.  —  An  “  A. 


and 


-  y„w  -  viw 


BLOCK  V-4. 


Is  -ZR^  v1  s  2R? 


This  block  determines  whether  or  not  the  target  is  within  permissible  bounds: 
on  the  yw  axis.  R  is  the  burst  radius  of  the  cloud.  If  -2R  s  y*  s  2R,  then 
the  target  may  be  subject  to  the  vapor  cloud  effects.  If  y'  is  less  than  -2R 
or  if  y’  is  greater  than  +2R,  the  vapor  cloud  has  no  effect  upon  the  target. 


BLOCK  V-5.  Is  x1  >  -R? 


If is  greater  than  -R,  then  the  vapor  cloud  may  affect  the  target.  If 
x'  is  not  greater  than  -R  (that  is,  x’  is  less  than  -R),  the  vapor  cloud  has 
no  effect  upon  the  target. 
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COORDINATE  SYSTEM 
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BLOCK  V-6. 
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Is  1  "iMAX? 


If  i  is  equal  to  ij^^,  the  effects  of  the  vapor  cloud  for  eacli  of  the  impacts 
have  been  considered  upon  a  specific  target.  If  i  is  not  equal  to 
(that  is,  i  is  less  than  i^AX^'  then  the  effects  of  additional  impacts  must 
be  considered  for  the  target  under  consideration. 

BLOCK  V-7 .  Increase  i  by  One 

The  value  of  i  is  increased  by  one,  so  that  the  vapor  effects  of  the  next 
impact  upon  the  target  element  may  be  evaluated. 

BLOCK  V-8.  Compute  Time  Dependent  Vapor  Casualties  for  a 

Target  Element  and  Print  Out 

The  effects  of  the  various  dosages  of  vapor,  which  are  obtained  by 
breathing  and  gasping,  are  summarized  for  each  man  and  printed  out. 

Also,  there  is  a  print-out  for  each  man  which  gives  the  number  of  times, 
if  any,  in  which  he  received  a  casualty-producing  dose. 

BLOCK  V-9.  Isn  =  nMAX? 


If  n  is  equal  to  nj^^,  the  effects  of  each  of  the  vapor  clouds  for  all 
impacts  upon  each  of  the  target  elements  have  been  evaluated  and  pertinent 
information  printed  out;  therefore,  the  computer  enters  the  routine  for 
evaluating  the  effects  of  the  time  dependent  impingement.  If  n  is  not  equal 
to  njyj^x  ^hat  is,  n  is  less  than  njyjj^),  the  next  target  element  is  considered. 

BLOCK  V-10,  Increase  n  by  One;  set  1  =  .l 

The  value  of  n  is  increased  by  one,  so  that  the  vapor  cloud  effects  upon 
the  next  target  element  will  be  evaluated  for  each  of  the  impacts. 

BLOCK  V- 1 1 .  Have  Breathing  Parameters  Been  Generated? 

If  the  breathing  parameters  have  been  generated,  the  computer  enters  the 
dose  subroutine.  If  the  breathing  parameters  have  not  been  generated,  the 
computer  enters  the  breathing  parameter  subroutine. 
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Breathing  Parameters 


BLOCK  B-l.  Select  at  Random  tw. 


The  random  variable,  is  the  time  relative  to  the  beginning  >f  a  man's 
breathing  cycle  ar  impact  time.  The  variable  t, w  is  expressed  as  an  inte¬ 
gral  number  of  intervals  of  width  6.  The  inhale  time  is  2a 6;  the  exhale 
time  is  2b6;  a  and  b  are  integers. 

BLOCK  F-2.  Select  at  Random  Response  Delay  Time 

This  is  the  time  which  the  target  element  takes  before  reacting  to  the 
impact.  This  delay  time  is  expressed  as  an  integral  number  of  intervals 
of  width  6. 

BLOCK  B-3.  Store  Response  Delay  Time 

The  value  of  the  response  delay  time,  t  , ,  is  stored.  Note:  In  the  de¬ 
termination  of  the  breathing  parameters,  it  is  convenient  to  have  two  time 
scales.  One  time  scale  is  measured  from  impact  time.  The. other  scale 
is  measured  from  the  beginning  of  the  breathing  cycle  of  a  particular  target 
element.  This  time  scale  is  denoted  by  time  symbols  containing  a  prime 
(that  is,  t').  A  time  scale  in  terms  of  the  impact  is  shifted  ty,  time  units 
from  the  other  time  scale. 

BLOCK  B-4 .  Determine  t  , ' 


Symbol  t.^'  is  used  to  denote  the  time  relative  to  the  beginning  of  the 
breathing  cycle  when  the  target  element  begins  to  respond  to  the  impact. 
This  time  value  is  obtained  as  follows: 

t-l’  =  t-i  +  V 

BLOCK  B-5.  Select  at  Random  Masking  Times 

The  time  required  for  each  target  element  to  mask,  £tm,  is  determined  from 
a  distribution  of  masking  times,  which  are  expressed  as  an  integral  number 
of  inter  valsofwidth  6. 


BLOCK' Bt6.  Compute  and  Store  Masking  Times 

The  times  at  whichmaskingoccursis  determined  for  both  time  scales  from 
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the  following  equations: 

Sn  “  Vi  +  and 
tm’  =  +  tw- 

Breathhold  Rules 

There  are  five  blocks  which  are  used  to  determine  the  time  at  which  the 
target  element  begins  the  breathhold  procedure.  See  Figure  5.  The  rules 
which  are  used  in  establishing  the  time  for  breathhold  are  as  follows: 

(1)  If  at  the  end  of  the  response  delay  time,  a  man  is  in  the  inhale  portion 
of  the  cycle ,  he  will  continue  inhal  ig  until  he  reaches  the  cycle  peak , 


(2)  If  at  the  end  of  the  response  delay  time,  the  man  is  in  the  exhale 
portion  and  has  exhaled  less  than  50%  of  the  volume  of  air,  he  will  begin 
the  breathhold  immediately. 

(3)  If  at  the  end  of  the  response  delay  time,  the  man  is  in  the  exhale 
portion  of  the  cycle  and  has  exhaled  more  than  50%  of  the  volume  of.  air, 
he  will  continue  exhaling  and  inhale  50%  of  the  normal  volume  before 
breathhold  begins . 


BLOCK  B-7. 


Is  t  s  t  '? 
-1.  P 


If  t_j'  is  equal  or  less  than  tp'  then  the  man  was  in  the  inhale  portion 
of  the  cycle.  Otherwise,  the  man  was  not  in  the  inhale  portion  of  the 
first  cycle. 


BLOCK  B-8. 


set:tHi'  ~y 


In  accordance  with  rule  one,  as  stated  above,  the  man  begins  breathhold 
at  the  peak  of  the  cycle,  if  he  is  in  the  inhale  portion  at  the  end  of  the 
response  delay  time. 

t  '  +  t,  * 


BLOCK  B-9. 


is  y 5 


t  '  .  t,» 

If  t  'is  equal  to  or  less  than  _ ,  then  the  man  was  in  the  exhale 

“  *■  w 

portion  of  the  cycle  and  less  than  50%  of  all  the  volume  of  intake  air  had 
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V  +  li 

been  exhaled.  If  t  '  is  not  equal  to  or  less  than  r'  ■  ^ — ,  the  man  Is 
not  in  the  first  half  of  the  exhale  portion  of  the  cycle  ir.  'Which  the  impact 
occurred. 

BLOCK  B-10.  Set  t  '  =  t  .  ’ 

ill  - 1 


In  accordance  with  the  breathhold  rules ,  when  the  man  is  in  the  exhale 
portion  of  the  cycle  and  less  than  50%  volume  has  been  exhaled,  he 
begins  breathhold  immediately. 

BLOCK  B-ll.  Is  L  tj'? 


If  t  '  is  equal  to  or  less  than  t^\  then  the  man  was  in  the  exhale  portion 
of  the  cycle  and  more  than  50%  of  volume  had  been  exhaled.  If  t  ,  ’  is 
not  equal  to  or  less  than  t^',  then  the  man  was  not  in  the  breathing  cycle 
in  which  the  Impact  occurred. 

BLOCK  B-12.  Set  t  '  =  t. '  + 

_ Til _ 1 


In  accordance  with  the  third  breathhold  rule,  when  the  man  is  in  the  exhale 
portion  of  the  cycle  and  more  than  50%  of  his  intake  volume  has  been  ex¬ 
haled,  the  time  at  which  he  begins  breathhold  occurs  after  he  has  completed 
the  exhale  and  half  of  the  next  inhale. 

BLOCK  B-13.  1st  ,'st  ,'  +  t  '? 

-lip 


If  t_^'  is  equal  to  or  less  than  t  '  +  t  ',  then  the  man  is  in  the  inhale 
portion  of  the  second  cycle.  If  t_ ^ *  h  not  equal  to  .or  less  than  t  ’,  then 
the  man  is  not  in  the  inhale  portion  of  the  second  cycle. 

BLOCK  B-14.  Set  t  '  -  t.'  +  t  ' 

Ill  1  P 


In  accordance  with  the  first  breathhold  rule,  the  time  at  which  breathhold 
begins  is  at  the  peak  of  the  inhale  cycle. 

3ti*  +  tD' 

BLOCK  B-15.  Is  t.j'  *  ? 

3ti‘  +  tp« 

If  t,  is  equal  to  or  less  than - '  •  t*ie  r«‘isn  was  in  the  exhale  portion 
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of  the  second  cycle  and  less  than  50%  of  volume  had  been  exhaled.  If 

t  is  not  equal  to  or  less  than  3ti  +  *p  ,  then  the  man  is  not  in  the  first 
-1  - 2 - 

half  of  the  exhale  portion  of  the  second  cycle  since  the  impact  occurred. 
BLOCK  B-16.  SettHl'=t-l' 


In  accordance  with  rule  three,  when  the  man  is  in  the  exhale  portion  of 
the  cycle  and  less  than  50%  of  volume  has  been  exhaled,  the  breathhold 
is  begun  immediately. 

y 

BLOCK  B- 1 7 .  Set  tH1 '  =  2tx '  + 


In  the  logic  utilized  in  the  program,  it  is  implicit  that  the  man  will  be  in 
either  the  first  or  second  breathing  cycle  when  he  begins  to  respond  to 
the  impact.  Consequently,  with  a  negative  response  to  Block  B-15,  the 
man  is  in  the  exhale  portion  of  the  second  cycle  and  more  than  50%  of  the 
volume  of  air  has  been  exhaled.  Therefore,  in  accordance  with  the  third 
breathhold  rule,  he  will  complete  the  exhale  procedure  and  inhale  50%  of 
the  normal  inhale  volume  prior  to  holding  his  breath. 

BLOCK  B- 18.  Set  J  =  1 

The  following  blocks  are  concerned  with  determining  the  times  at  which 
gasps  following  the  first  breathhold  time  occur. 

BLOCK  B- 1 9 .  Select  at  Random  a  Breathholdinq  Time 

From  a  probability  distribution  of  breathholding  times  a  value  denoted  as 
BHTj  is  obtained.  This  variable,  BHTj,  is  expressed  as  an  integral  number 
of  intervals  of  width  6. 

BLOCK  B-20.  Store  tQ.' 


The  time  at  which,  for  example,  the  first  gasp  occurs,  is  equal  to  the  time 
of  the  first  breathhold  plus  the  breathhold  time.  (See  Figure  5.)  Note  that 
the  breathhold  time  includes  the  time  for  exhaling. 

BLOCK  B-21.  Is  tr  '  <  t  '? 

Gj  m 


This  block  inquires  whether  or  not  the  j**1  gasp  occurred  prior  to  the  time 
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at  which  masking  was  completed.  If  tQj'  is  less  than  tm\  then  the  man 
has  not  masked  himself  at  the  beginning^  of  the  jth  gasp.  If  tGj'  is  not 
less  than  t  ',  masking  has  occuired  prior  to  the  beginning  of  the  jtl1  gasp; 
consequently,  the  jtfl  gasp  does  not  produce  any  vapor  dose. 


BLOCK  B-22 . 


Increase  i  by  One 


The  test  of  Block  B-21  is  applied  to  the  (j  +  l)th  gasp. 


BLOCK  B-23. 


Select  at  Random  BHT, 


From  a  probability  distribution  of  breathhold  times ,  a  value  for  the  next 
breathhold  time  will  be  selected. 


BLOCK  B-24 . 


Compute  t„.' 

Ca) 


All  values  of  t  '  for  j  equal  to  2, 3, ... ,  are  obtained  by  adding  to  the 
time  at  which  the  previous  gasp  began,  the  gasp  duration,  GD,  and  the 
next  breathholding  time . 


BLOCK  B-25, 


Set  tGj '  =  oo 


When  it  is  determined  that  the  masking  time  occurred  prior  to  the  gasping 
time  a  large  number,  9.9  x  1040,  is  placed  in  storage  to  represent  infinity. 


BLOCK  B-26. 


Select  at  Random  TV 


The  volume  of  air  intake  is  selected  at  random  from  a  probability  distribution 
of  tidal  volumes .  Note;  The  following  breathing  parameters  have  been- 
generated  for  a  specific  target  element: 

W 

‘gi’-  tcz’"**'00 


BLOCK  B-27 . 


Exit  to  Dose  Subroutine 


After  the  breathing  parameters  have  been  determined,  the  computer  begins 
the  dose  subroutine. 


Dose  Routine 
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BLOCK  R-0.  Compute  tA  and  tD 


This  block  computes  the  times  of  arrival  and  departure  of  the  gas  cloud 
for  a  particular  target  from  the  following  equations: 


and 


x!  +  2R 
D~  Q 


If  the  formula  leads  to  negative  tA,  tA  is  set  equal  to  zero.  Both  t^  and 
tD  are  expressed  as  an  integral  number  of  intervals  of  width  6. 

BLOCK  R-l .  Compute  t^'  =  t^  +  tyy 


This  block  computes  the  time  of  arrival  of  the  gas  cloud  with  respect  to 
the  beginning  of  the  man's  breathing  cycle. 

BLOCK  R-2.  Is  tj' «t  tm'  ? 


If  the  time  at  which  the  gas  cloud  arrives  is  greater  than  the  time  at  which 
the  man  puts  on  his  gas  mask,  then  the  man  receives  no  dose.  The  effects 
of  the  next  impact  will  then  be  determined  with  respect  to  this  same  man. 

BLOCK  R-3.  IstE'atHl’? 


If  the  time  at  which  the  gas  cloud  arrives  is  equal  to  or  greater  than  the 
time  of  first  breathhold,  then  the  machine  will  determine  the  gasp  dose. 
However,  if  the  time  at  which  the  gas  cloud  arrives  is  less  than  the  time 
of  first  breathhold  the  steady-breathing  dose  will  be  determined. 

BLOCK  R-4.  t  '-t* 

Compute  ~  —  :=ymx 


This  block  determines  the  number  of  increments  which  will  be  used  in  the 
numerical  integration  for  ascertaining. the  quantity  of  vapor  inhaled 
prior  to  the  beginning  of  the  first  breathhold. 


BLOCK  R-5. 


Set  y  -  1 
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The  index,  y,  is  used  to  denote  the  number  of  the  next  numerical  integration 
to  be  computed.  The  time,  t,  is  set  at  the  midpoint  of  the  first  interval 
after  the  gas  cloud  arrives . 


BLOCK  R-6, 


Determine  q  =  MOD 


The  breathing  cycle  increment  number  for  the  target,  corresponding  to  the 
time  at  which  the  gas  cloud  arrives,  is  determined. 

BLOCK  R-7.  Look-up  £f(q) 

The  fraction  of  TV  inhaled  in  the  q4*1  increment  of  the  breathing  cycle  is 
determined  from  a  table  placed  in  the  machine. 

BLOCK  R-8.  Is  4f(q)  >  Q? 

If  the  fraction  of  TV  inhaled  is  greater  than  zero,  the  incremental  volume 
inhaled  is  obtained  as  follows: 

*v(q)  =  TVUf(q)] 

If  the  quantity  of  vapor  inhaled  is  equal  to  zero,  this  computation  is 
avoided. 

BLOCK  R-9.  Is  t  =  tA  +  (Zy  -  1)  -|-  z  tD? 


If  the  mid-point  of  the  interval  is  greater  than  the  time  of  departure  of 
the  gas  cloud,  then  no  vapor  is  inhaled  since  the  gas  cloud  has  passed 
by.  If  the  mid-point  of  the  interval  is  less  than  the  time  at  which  the 
gas  cloud  departs ,  then  the  man  may  breathe  while  the  gas  cloud  is 
around  him. 

BLOCK  R-10.  Form  CAv(q)]  Lx(t,x',  y‘)]  and  add  to  Bn 


If  the  gas  cloud  has  not  passed  by,  then  the  quantity  of  vapor  inhaled  is 
computed  as  a  product  of  the  quantity  of  air  inhaled  times  the  concentration 
of  the  vapor.  The  concentration  of  the  vapor  is  a  function  of  the  time  since 
the  impact  and  the  x  and  y  coordinates  in  terms  of  the  wind-coordinate 
system.  This  product  is  then  added  to  the  previous  quantity  of  steady¬ 
breathing  dosage  for  the  man  which  is  stored  in  Bn* 


41 


BLOCK  R-ll. 


Isy=yMAX? 


If  y  is  equal  to  7j all  of  the  intervals  in  which  the  man  could; have 
breathed  the  vapor  prior  to  beginning  the  breathhold  routine  have  been 
considered.  The  machine  will  next  determine  the  quantity  of  vapor  inhaled 
during  the  gasping  routine.  If  y  is  not  equal  to  y^^x*  (that  is,  y  is  less 
than  y^p^)i  then  the  subsequent  interval  will  be  considered, 

BLOCKR-12.  Increase  v  by  one;  (t  by  61 

7  and  t  are  incremented  by  one  and  6,  respectively,  for  use  in  determining 
the  concentration  of  gas  in  the  next  interval. 

BLOCK  R-13.  Ist>tm? 


If  t  is  greater  than  tm,  then  the  man  is  masked  prior  to  the  mid-point  of 
the  subsequent  interval  and,  as  a  result,  he  will  receive  no  vapor.  The 
machine  then  will  determine  the  effects  of  the  next  Impact  upon  this  man. 

BLOCK  R-14.  Is  q  =  q^^X?  (Q^AX  ~  "I* ) 

If  q  Is  equal  to  q^^,  then  the  dose  for  ali  intervals  in  this  breathing 
cycle  has  been  computed  for  this  man.  If  q  is  less  than  q^^,  then  the 
man  is  still  in  the  same  breathing  cycle. 

BLOCK  R- 1 5 .  Increase  q  bv  One 

The  index  q  is  increased  by  one  so  that  the  fraction  of  the  tidal  volume 
Inhaled  during  the  subsequent  interval  can  be  obtained  as  in  Block  R-7. 

BLOCK  R" i  6 .  Set  g  =  1 

The  index  q  is  set  equal  to  one  so  that  the  next  breathing  cycle  for  the 
man  can  be  properly  initiated. 

BLOCK  R- 17.  Set  \  =  1 

The  index  for  the  gasp  time  is  set  equal  to  one. 

BLOCK  R-18.  Is  tGj*  =  co? 
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*  ^*"w*»«  -a^*— — ! u  4riea*i*'««miii»f*i»i*Maiii'tf^^ 


If  t  ‘  is  equal  to  oo,  this  implies  that  the  man  had  his  gas  mask  on  when 
he  m^de  the  j**1  gasp.  If  w  is  not  equal  to  oo,  the  man  did  not  have  his 
gas  mask  on  and  the  quantity  Oi  vapor  inhaled  will  be  computed  later. 


BLOCK  R-19. 


Set  y  =  1 


The  interval  index  y  (not  the  same  y  as  previously  used)  is  set  equal  to 


BLOCK  R-20. 


Form  tGj=tGj*"tw 
Set  t  =  L, 


The  time  for  the  j*-*1  gasp  is  obtained  in  terms  of  time  with  respect  to  the 
impact  time. 


BLOCK  R-21 . 


Is  tGj  +  GD  £  tA? 


If  tGj  +  GD  is  equal  to  or  less  than  the  time  of  arrival  of  the  vapor  cloud, 
then  no  vapor  has  been  inhaled  during  the  jth  gasp;  otherwise,  during  the 
jtn  gasp  some  of  the  gas  may  have  been  inhaled  by  the  man. 


BLOCK  R-22. 


Increase  j  by  One 


The  index  j  is  increased  by  one,  so  that  the  effects  during  the  subsequent 
gasp  can  be  determined. 


BLOCK  R-23. 


tGj  + (2y  -  1)  Y  *  *D? 


If  the  time  of  the  mid-point  of  the  y**1  interval  during  the  jt  .gasp  is  equal 
to  or  greater  than  the  time  of  departure  of  the  gas  cloud,  then  no  vapor 
was  inhaled  during  the  y**1  interval,  otherwise  some  vapor  may  have  been 
inhaled  during  the  yth  interval. 


BLOCK  R-24. 


IS  tGj  +  (2y-l)f  >tm? 


If  the  time  of  the  mid-point  of  the  yth  interval  on  the  jtn  gasp  is  greater 
than  the  time  of  masking,  then  no  vapor,  was  inhaled  during  the  y**1  interval. 
The  effects  of  the  next  impact  are  then  ascertained,  otheiwise,  some  vapor 
may  have  been  inhaled  during  the  yt^1  interval. 


BLOCK  R~25, 


Form  Avgy  X(tGj  +  (2y  -  1)  ~)  and  add  to  (Gn) 


Vapor  dose  in  the  y1*1  interval  of  the  j**1  gasp  is  obtained  by  obtaining  the 
product  of  the  volume  inhaled  and  the  vapor  concentration. 

BLOCK  R-26.  Is  y  =  ? 


If  y  is  equal  to  GD  divided  by  6,  then  the  complete  vapor  dose  due  to 
the  jth  gasp  has  been  ascertained  and  the  subsequent  gasp  will  be  con¬ 
sidered. 

BLOCK  R- 27.  Increase  y  by  One ,  t  by  8 

The  index  y  and  time  t  are  incremented  by  one  and  d.respectively,  in  order 
to  evaluate  the  vapor  dose  during  the  next  interval. 
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|  IV.  COMPUTER  PROGRAM 

I 

4.1  This  section  contains  the  literal  computer  programs  in  ALGOL  • 
and  FORTRAN.  These  are  the  specific  forms  in  which  the  simulation  is 
read  into  the  computer. 

4.2  At  present  a  detailed  annotation  of  the  programs  relating,  them  to 
the  blocks  in  the  flow  charts  has  not  been  made.  This  will  be  carried  out 
in  the  next  contract  year. 


SIMULATION  MODEL  OF  THE 
1 55-mm  HOWITZER  WEAPON  SYSTEM 


THE  ALGOL  VERSION  OF  THE  COMPUTER  SIMULATION 
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